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1. INTRODUCTION 

The economic dispatch (ED) is one of the optimization problems in power system operation and 
planning to allocate the shared power demand between the generating units. Thus, optimal power system 
operation is important in electrical networks to ensure the system can operate at minimal cost. Therefore, the 
aim of the ED is to minimize the total cost of generation and satisfy the system and operating constraints [1]. 

Various optimizations methods have been proposed and applied to solve ED problem over the latest 
decades and can be classified into two main categories such as classical method and heuristic method [2]. 
The classical methods such as newton's method [3], quadratic programming technique [4], interior point [5], 
lambda iteration method [6], evolutionary programming (EP) techniques [7] and dynamic programming [8] 
are widely used for solving convex and smooth cost function of ED problem. However, most of these 
methods have difficulty for solving nonconvex or nonsmooth problems. To solve this problem, 
the non-conventional or heuristic methods are developed to solve the complicated and highly nonconvex 
optimization problem. The optimization algorithms such as ant colony optimization (ACO) [9], artificial bee 
colony (ABC) [10], firefly algorithm (FA) [11], particle swarm optimization (PSO) [12], teaching—learning- 
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based optimization (TLBO) [13], genetic algorithm (GA) [14] and adaptive charged system search algorithm 
(ACSS) [15] have been solved the complex ED problem. The performance of these algorithms is better than 
classical optimization methods in many aspects for instance fast, robust and easy to adjust according to the 
problem. In some application, these algorithms are suffering from slow convergence rate, stuck at local 
solution and required proper parameter tuning to obtain optimal solution. 

Furthermore, the hybrid method has been introduced by combining two or more algorithms in 
order to mitigate their weaknesses and use their strengths to provide better performance for solving 
optimization problems [16]. The proposed hybrid algorithms such as GA-PS-SQP [17], NM-FAPSO [18] 
and differential evolution algorithm-PSO [19] shows highly efficient technique to solve the ED problem. 
However, it required long computational time and complex programming since two or more 
algorithms are used. 

Recently, a new nonconventional algorithm called the LSA algorithm [20] has never applied for 
solving ED problem. It was successfully used and implemented to optimize the different applications such as 
binary optimization [21], nuclear reactor controller [22], and Fuzzy Logic PV Inverter Controller [23]. 
It shows the LSA can provide superior performance over other algorithms. Therefore, this paper proposed 
LSA as a new approach for solving ED problems with VP. The performance of LSA has been compared with 
other common methods. The remaining of this paper is organized as follows: Section 2 provides the 
Mathematical formulation of ED problem considering VPE and losses, section 3 presents the LSA technique, 
Section 4 introduces the simulation results for the test system and its analysis. 


2. PROBLEM FORMULATION OF ECONOMIC DISPATCH (ED) 
2.1. Objective function 

The main aim of any ED problem is minimizing the system operation cost in order to fulfill the 
power demand alongside the generator limits. The fuel cost function can be differentiated as a simplified cost 
function and modified function by including with VPE. The cost function of production units can be 
represented by a quadratic function as follows: 


Ng 
F(Fe)= Fi (Fai) (1) 
F, (Po;) = @,Po, + O.Pey + 6; (2) 


Considering VPE for ED problem in the fuel cost functions exhibit a bigger variation in the 
generating units with multi-valve steam turbines. The generated power will be changed when every steam 
valve closes or opens. In the curves of the heat ratio introduces ripples by VPE. The cost function of ED will 
add a sinusoidal term and it can be defined mathematically as: 


F(P)=(a, +b, P,, +eP’g ) +e, sin( F (Penn ~ Pav) (3) 


where F(Pg) is the total cost of production, F(Pgi) is the unit function of fuel cost i; a;,b;,c;,e;, and F; are the 
unit coefficients of fuel cost 7; Pgi is the unit output i of real power. 


2.2. System operation constraints 

It should satisfy the following system constraints when ED Optimization carried out. The constraints 
considered for this research are discussed in sections 2.2.1 and 2.2.2. 
2.2.1. Equality constraint 

The equality constraint is represented the power balance equation and it is used in the optimization 


to ensure that the total power generated must be equal to the power demand in the case without losses in the 
system can be written as: 


y P= 6, (4) 
i=l 


If consider power losses, total generated power must be equal power losses and power demand as follows: 
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DP HPP, (5) 


i=l 


where P; is the power generated, P; Indicates transmission loss, Pp is the load demand and n is the overall 
amount of the generating units. The power losses can be calculated as: 


lag SP Bi Pa + 2 Bu, vi + Boo (6) 


i=l i=l 
where Bj, Bi, and Boo are to determine the loss coefficient matrix 


2.2.2. Inequality constraint 
A power limit constraint is used to ensure that the generating unit operates within the minimum and 
maximum limits. This constraint presented as follow: 


P.,, <P. <P (7) 


where Pinin and Pinax are minimum and maximum generating power limits of i” generating units. 


3. LSA ALGORITHM 

The LSA algorithm is a new metaheuristic method developed in 2015 by [20]. This method can be 
used to optimize complex nonlinear problems. It derived from a natural lightning phenomenon as a step 
leader propagation mechanism [24]. 

Nearby the thundercloud region can be found oxygen and nitrogen and hydrogen molecules. 
During the water molecules are freezing within a thundercloud, parts of water molecules are unable to fit the 
ice structure. Thus, these molecules at high speeds will be ejected from the forming ice. Therefore, 
the hydrogen and oxygen atoms are detached and ejected randomly in different directions as projectiles. 
These projectiles travel over the atmosphere and start the ionization path over collision and transition into the 
step leader. In this proposed algorithm, every projectile will create a step leader or channel that represents the 
initial population size. The projectile concept in this algorithm is highly similar to “particle” used in PSO. 

The projectiles are taken into consideration which is known as the fast particles involvement is 
recognized in the structure formation of the binary tree of the step leader. Also in the formation of concurrent 
of two leader tips at points of the fork instead of the conventional mechanism of step leader which utilizes the 
streamers concept [20]. There are three types of projectiles that are developed such as the transition 
projectiles to produce initial solutions, space projectile provides exploration for good solutions and the lead 
projectile performs exploitation. Figure | presents a flowchart of proposed LSA for solving ED problems. 


D 
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Figure 1. LSA Flowchart of proposed LSA for implementation on ED problem 
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3.1.1. Transition projectiles 

This projectile type is responsible for producing initial step leaders form the early stage because the 
transition projectile ejected from thunder in direction randomly. Therefore, by using a uniform probability 
distribution, this stage can be modeled in a random number. It can be described as follows: 


- fora<x’ <b 
f(x")=4 ba (8) 


0 forx<aorx’ >b 


where x" represents the random number trying to give the solution or provide initial tip energy E*/ of step 
leader S'i and the population of N can be given as S* = [s!1, s!2, s!3, . ., sIN], and P™ = [p™!, p7?, p™?,.., p7] 
for satisfying the dimension of the required solution. 


3.1.2. Space projectile 

Space projectile tries to explore the process in order to achieve the best position of leaders through 
ionizing the region of the old leaders. Therefore, the step leaders will be changed after the ionization of step 
leader tips tip in stept+1. The space projectile P* location = [p*!, pS’, p®?, . , pSX] can be slightly designed as a 
random number produced from an exponential distribution with shaping parameter yu. The following equation 
can be expressed as: 


. 4 tis for x° >0 
f(x )=) a (9) 
0 


for x° <0 


where x5 presents a random variable. The location or direction of space projectile can be clearly denoted in 
the equation and can be in charge of by shaping parameter 1. The distance between space projectile p* and 
lead projectile p“ is taken by pi for specific p*. According to this description, p* position at step+1 
is given as: 


Pe, =P £exprand(u,) (10) 


where exprand is an exponential variable created randomly. When p“ has a negative value, then it should 
subtract the produced random number due to (10) that provides only positive values. Therefore, 
new position p* yew does not ensure stepped leader propagation until it can find a good solution. 


3.1.3. Lead projectile 

The lead projectile travels nearest to the ground and its associated projectile does not have enough 
potential for the large ionizing section in front of the leading tip. Therefore, it can obtain from normal 
distribution a modeled the lead projectile as a random number. 


f (x") = 1 eo -p)" 120° (11) 


~ oN2a 


where pL. is a shaping parameter to specify the present location of lead projectile p’ and o is the scaling 
parameter that is responsible to exploit the projectile ability. Thus, p', o, and ps are decreasing exponentially 
and progress to discover the best solution. Based on this idea, at step+1 the p" position can be written as: 


Pe. =p" +normrand (1, ,0,) (12) 
where normrand presents a random number produced by the distribution function. p* is an updated location 
of the lead projectile while p’ new is the new location. In a similar way, if p’ new does not get a good 
solution, then step leader will provide to a new location. 


4. SIMULATION RESULTS 
In this research, the LSA algorithm is applied to two different test systems as discussed in section 
4.1.1-4.1.2, to validate the proposed algorithm for ED problems. In order to investigate the performance of 
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the LSA algorithm, the achieved results have compared with different experimental results of the latest 
literature. The power balance and generator limits constraints are considered for all test systems. 
LSA technique has implemented in Matlab R2018b considering 6 and 13 units from the previous studies. 
The parameters used in the experimental study of LSA; 50 projectiles as the population size of the algorithm, 
the maximum number of the iteration is 500 and 50 multiple runs in every case study having 50 various 
initial solutions. 


4.1.1. Test case 1 

This test system has implemented for six thermal generating units using the proposed 
LSA approach. The demand for power in this test case is 283.4 MW and the transmission loss coefficient is 
given in [25]. The cost function with the VPE and losses are considered. Table | presents the comparisons of 
optimal results obtained by the LSA with existing approaches i.e. MSG-HS [26], PSO [26], and GA [25]. 
The total value of generation cost that can be achieved by LSA reached 925.122 ($/h) is a bit lower and better 
than compared methods. It provides also other statistics results for LSA performance i.e. simulation time, 
best, average and the worst cost. The convergence and robustness of the LSA for minimizing cost is shown in 
Figure 2(a-b). Figure 2(a) it is observable that the LSA is trying to converge the cost function with the best 
iteration from 40-500" for cost value. Moreover, Figure (b) illustrates the consistency of the optimal result 
after 50 runs with standard division around 12.02. It found that 31 runs can produce minimum cost results of 
about 925.122 $/hr. 


Table 1. Comparison results for 6 units 


Units Power output (MW) 
LSA MSG-HS [26] PSO [26] GA [25] 
Pl 199.628 199.633 197.865 150.724 
P2 20.074 20 50.337 60.870 
P3 24.239 23.762 15 30.896 
P4 20.073 18.393 10 14.214 
PS 17.912 17.102 10 19.489 
P6 12.506 15.692 10 15.915 
The power generated (MW) 294.37 924.583 295.758 292.11 
Losses (MW) 10.97 11.183 11.802 8.706 
Best Cost ($/hr) 925.122 925.640 925.758 996.037 
Average Cost ($/hr) 930.88 925.64 925.76 - 
Worst Cost ($/hr) 962.959 928.6 928.6 1117.13 
Power Demand 283.4 MW 283.4 MW 283.4 MW 283.4 MW 
Simulation time (s) 1.285 0.621 0.353 0.578 
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(a) Convergence characteristic 


(b) Robustness of Cost 


Figure 2. LSA performance for test case | 


4.1.2. Test case 2 

A system of 13 generators considering VPE has examined with two different load demands i.e. 1800 
MW without transmission losses and 2520 MW with transmission losses. The input data for test case 2 with 
and without losses have taken from [1]. The obtained result from 1800 MW load demand have compared 
with other algorithms as improved coordinated aggregation-based PSO (ICA PSO) [27], quantum-inspired 
particle swarm optimization (QPSO) [28], chaotic differential evolution and sequential quadratic 
programming (DEC SQP) [29] as shown in Table 2. While, the test results of 2520 MW load demand 
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considering losses are also compared with existing solutions i.e. OIWO [1], SDE [30], ORCCRO [31] as 
shown in continuous para of Table 2. It is observable that the minimum cost obtained using LSA is about 
17963.95 ($/hr) for without loss test and the lowest cost of about 24507.66 ($/hr) for losses test. It also 
provides other statistics results i.e. simulation time, best, average and the worst cost for LSA in comparison 
with losses test system literature. 

The convergence and robustness of the LSA for minimizing cost is shown in Figure 3(a-d). 
Figure 3(a) represents the convergence characteristics of the 1800 MW load demand using LSA. The cost 
value has converged for the cost function with the best iteration from 150-500". Moreover, Figure 3(b) 
illustrates the consistency of the optimal result after 50 runs with a standard division of 50.86. Similarly, 
Figure 3(c) for 2520 MW load demand the converged value tor cost function with the best iteration is 180%- 


500" and the robustness of load demand with losses is shown in Figure 3(d). 


Table 2. Results comparison for test case 2 


Units Without losses (Pp=1800 MW) With losses (Pp=2520 MW) 
LSA ICA PSO [27] QPSO [28] LSA OIWO [1] SDE [30] | ORCCRO [31] 
Pl 628.31 628.32 538.56 628.3185 628.3185 628.32 628.32 
P2 224.16 149.96 224.70 299.1993 299.1989 299.20 299.20 
P3 148.18 299.75 150.09 295.1809 299.1991 299.20 299.20 
P4 109.86 109.86 109.87 159.7331 159.7331 159.73 159.73 
P5 109.86 109.86 109.87 159.7331 159.7331 159.73 159.73 
P6 109.86 60 109.87 159.7331 159.7331 159.73 159.73 
P7 60 109.87 109.87 159.7331 159.7330 159.73 159.73 
P8 109.86 109.87 109.87 159.7331 159.7331 159.73 159.73 
P9 109.86 109.87 109.87 159.7331 159.7330 159.73 159.73 
P10 40 40 T7 Al 77.39991 77.3953 77.40 77.40 
Pll 40 40 40 114.7998 113.1079 113.12 113.14 
P12 55 25 55.01 92.39991 92.3594 92.40 92.40 
P13 55 55 55.01 92.39991 92.3911 92.40 92.40 
Power generated(MW) 1800 1800 1800 2558.097 2560.368 2560.44 2559.44 
Losses - - - 38.097 40.3686 40.43 39.43 
Best cost ($/hr) 17963.95 17964.06 17969.02 24507.66 24514.83 24514.95 24513.99 
Worst cost ($/h) 18167.54 - - 24692.31 24514.83 - 24513.91 
Average cost 18078.49 - - 24556.99 24514.83 24516.31 24513.91 
Solution time (s) 9.4 - - 14.3 5.3 - 8 
- not reported by the authors 
18700 fy 18200 —> . r 7 
18600 + 
jasoel: 18150 + 4 
a 18400 + S| € 18100 4 
18200 + Py 
18100 + 18000 - se | 
18000 + 4 ei | 
eee vr T00 200 300 100 $00 ; : , ; 
Number of Iterations ° ia Ronee of eee . a 
(a) (b) 
‘ 0° Number of Iterations ° ” NantecstRuns: 
(c) 


(dy 


Figure 3. LSA performance for test case 2. (a) Convergence characteristic without losses (b) Robustness of 
Cost without losses (c) Convergence characteristic with losses (d) Robustness of Cost with losses 
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5. CONCLUSION 

Minimizing the generation cost and profit increment is a priority in the planning and operation of 
power systems. Therefore, this paper proposed LSA as a new metaheuristic algorithm for solving the ED 
problems considering the nonsmooth cost function (valve point effect) and transmission loss. The proposed 
LSA has been tested on two test cases varied in complexity and size. The optimal results obtained by LSA 
has been compared with the results obtained by other methods that reported in the literature for the standard 6 
and 13 units test system. It found that LSA can produced minimum cost of 925.122 $/h and 17963.95 $/h for 
6 units and 13 units system respectively. The convergence characteristic, robustness, simulation time of 
proposed LSA are also investigated in this paper for both test systems. The simulation results proved that 
proposed LSA is outperformed than compared method in term of generation cost. Further investigation can 
be performed to evaluate LSA performances for minimizing cost and emission level of power generation and 
other optimization problems. 
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